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Abstract 

A method for finding exact solutions of nonfinear differential equations 
is presented. Our method is based on the application of the Newton 
polygons corresponding to nonlinear differential equations. It allows one 
to express exact solutions of the equation studied through solutions of 
another equation using properties of the basic equation itself. The ideas 
of power geometry are used and developed. Our approach has a pictorial 
rendition, which is is illustrative and effective. The method can be also 
applied for finding transformations between solutions of the differential 
equations. To demonstrate the method application exact solutions of 
several equations are found. These equations are: the Korteveg - de Vries 
- Burgers equation, the generalized Kuramoto - Sivashinsky equation, the 
fourth - order nonlinear evolution equation, the fifth - order Korteveg - 
de Vries equation , the modified Korteveg - de Vries equation of the 
fifth order and nonlinear evolution equation of the sixth order for the 
turbulence description. Some new exact solutions of nonlinear evolution 
equations are given. 

Keywords: power geometry, exact solution, nonlinear differential equation, 
traveling wave. 

PACS: 02.30.Hq - Ordinary differential equations 

1 Introduction 



One of the most important problems of nonlinear models analysis is the construc- 
tion of their partial solutions. Nowadays this problem is intensively studied. We 
know that the inverse scattering transform [1-3] and the Hirota method [3-5] are 
very useful in looking for the solutions of exactly solvable nonlinear equations. 
However most of nonlinear equations are nonintegrable ones. For constructing 
solutions of such equations the following methods are used: the singular mani- 
fold method [7-13], the Weierstrass function method [13-15], the tanh-function 
method [16-21], the Jacobian elliptic function method [22-24] and the trigono- 
metric function method [25,26]. Most of these methods proceed from an a priory 
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expression for an unknown solution. Usually this expression is a polynomial in 

elementary or special functions (for example, in trigonometric, Weierstrass, Ja- 
cobian and some other). As a result all solutions outside the given family are 
lost. This disadvantage significantly decreases the strength of such methods. 

Lately it was made an attempt to generalize most of these methods and as 
a result the simplest equation method appeared [27, 28] . Two ideas lay in the 
basis of this method. The first one was to use an equation of lesser order with 
known general solution for finding exact solutions. The second one was to take 
into account possible movable singularities of the original equation. Virtually 
both ideas existed though not evidently in some methods suggested earlier. 

However the method introduced in [27, 28] has one essential disadvantage 
concerning with an indeterminacy of the simplest equation choice. This disad- 
vantage considerably decreases the eff'ectiveness of the method. In this paper 
we present a new approach for finding exact solutions of nonlinear differential 
equations, which greatly expands the method [27, 28] and which is free from 
disadvantage mentioned above. The method does not postulate the simplest 
equation but it allows one to find its structure using the properties of the equa- 
tion studied itself. Our method develops the ideas of power geometry [29-32] . 
With a help of the power geometry we show that the search of the simplest 
equation becomes illustrative and effective. Also it is important to mention 
that the results obtained are sufficiently general and can be applied not only to 
finding exact solutions but also to constructing transformations for nonlinear 
differential equations. 

The outline of this paper is as follows. The application of the Newton poly- 
gons of nonlinear differential equations and the method algorithm for finding 
exact solutions is presented in section 2. Our method is applied to look for 
exact solutions of the Korteveg - de Vries - Burgers equation, nonlinear evo- 
lution equation of the fifth order and the generalized Kuramoto - sivashinsky 
equations in sections 3, 4 and 5 accordingly. Self similar solutions of the fifth 
order Korteveg - de Vries equations are considered in sections 6 and 7. Solitary 
and periodical waves of sixth order nonlinear evolution equation for description 
of the turbulent processes are found in sections 8 and 9. 



2 Method applied 

Let us assume that we look for exact solutions of the following nonlinear n - 

order ODE 

Mn[y{z),y,{z),y,^{z),...,z]=0. (2.1) 

In power geometry any differential equation is regarded as a sum of ordinary 
and differential monomials. Every monomial can be associated with a point on 
the plane according to the following rules 

Ci z'^^y'^^ ^ (gi, 52), C2^^ {-k, 1). (2.2) 

Here Ci, C2 are arbitrary constants. When monomials are multiplied their 
coordinates are added. The set of points corresponding to all monomials of a 
differential equation form its carrier. Having connected the points of the carrier 
into the convex figure we obtain a convex polygon called the Newton polygon 



2 



of differential equation. Thus nonlinear ODE 1)2. 1|) can be characterized by 
the Newton polygon Li on the plane. The periphery of the polygon consists of 
vertexes and edges. By {rj^"*} let us denote the vertexes and by {r^-^''} the edges. 
Most of edges and vertexes of the Newton polygon define power or non - power 
asymptotics and power expansions for the solutions of the equation [29-32]. 

Now let us assume that a solution y{z) of the studied equation can be ex- 
pressed through solutions Y{z) of another equation. The latter equation is 
called the simplest equation. Consequently we have a relation between y{z) and 
Y{z) 

y{z)^FiYiz),Y,iz),...,z). (2.3) 

The main problem is to find the simplest equation. Substitution H2.3|l into 
basic equation H2.1(l yields a transformed differential equation, which is in its 
turn characterized by the polygon L2. By L3 denote the Newton polygon cor- 
responding to the simplest equation. Analyzing the Newton polygon L2 we 
should construct the polygon L3. Let the edge F^^^ belong to L2 and the edge 
f^^^ belong to L3. It can be proved that if r^^-* and f have equal external 
normal vectors, then corresponding asymptotics coincide accurate to the nu- 
merical parameter. Under the term external vector we mean the vector that 
goes out of the Newton polygon. Consequently, suitable Newton polygon L3 
has all or certain part of edges parallel to those of L2. Besides that when the 
Newton polygon L3 moves along the plane his apexes should cover the support 
of the transformed differential equation. From this fact we see in particular that 
suitable L3 is of equal or lesser area than L2. Let us suppose that we have found 
such polygon L^. Then we can write out the simplest equation 

E.raiY{z),Y,{z),...,z) = 0. (2.4) 

It is important to mention that the choice of the simplest equation is not unique. 
If the following correlation 

M„(F(r, y„ . . . , z)) = REraiY, Y,,...,z) (2.5) 

(where i? is a differential operator) is true then it means that for any solution 
Y{z) of the simplest equation 1)2. 4|1 there exists a solution of (|2.1I) . Generally 
speaking, any nonlinear solvable differential equation can be the simplest equa- 
tion. The only requirement is the following: the order of the simplest equation 
should be lesser than the order of the transformed differential equation. But the 
most important simplest equations are those that have solutions without mov- 
able critical points. If the general solution (or a partial solution) of the simplest 
equation can be found, then we get explicit representation for a solution of the 
equation studied. Otherwise we have only the relation 1)2. 3|l between solutions 
of and 123. 

The most useful examples of the simplest equations are the following: the 
Riccati equation 

Y, + Y'^~a{z)Y-b{z)=^0, (2.6) 
the equation for the Jacobi elliptic functions 

Rl = ^AR^ + aR^ + hR^ + cR + d, (2.7) 
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and the equation for the Weiershtrass elhptic functions 



Rl^-2R^ + aR^ + bR + c. 



(2.8) 



Solutions of equations H2.6|l . (|2.7(l . and 1)2.8(1 do not have movable critical 
points. 

Now we are able to state our method. It is composed of six steps. 

The first step. Construction of the Newton polygon Li, which corresponds 
to the equation studied. 

The second step. Determination of the movable pole order for solutions of 
the equation studied and transformation of this equation using the expression 



The third step. Construction of the Newton polygon L2 corresponding to 
the transformed equation. 

The fourth step. Construction of the Newton polygon L3, which will char- 
acterize the simplest equation. This polygon should posses properties discussed 
above. 

The fifth step. Selection of the simplest equation with unknown parameters 
that generates the polygon L3. 

The sixth step. Determination of the undefined coefhcients, which are present 
in the transformation ^Z.?t\ and in the simplest equation. 

Remark 1. Very often suitable simplest equations can be found without 
making transformation ((2.3(1 . In this case ((2.3(1 is an identity substitution and 
we set y{z) = Y{z), L2 = Li. 

Remark 2. In some cases the transformation can be included into the sim- 
plest equation. Then again ((2.3(1 is an identity substitution. 

Remark 3. The most powerful transformations, i.e. the transformations that 
generate new classes of exact solutions, are those that change the pole order of 



3 Exact solutions of the Korteveg — de Vries — 
Burgers equation 

To demonstrate our method application let us find exact solutions of the Ko- 
rteveg - de Vries - Burgers equation 



(E31). 



y{z). 



Ut + UUx + f3u. 



'XXX 



0. 



(3.1) 



It has the travelling wave reduction 



u(a;, t) 



w{z), 



z — X — C'ot, 



(3.2) 



where w{z) satisfies the equation 



E[W] = 



1_ Q 

— lywz + -w — Cqw + Ci = 0. 



(3.3) 



Here Co, Ci, a, and /3 are constants. In the case. 



Co 




(3.4) 
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equation H3.3|) is integrable and was solved by Painleve [5,6]. Its general solu- 
tion is expressed via the Jacobi elliptic function. However equation H3.3|l also 
describes a solitary wave. The solution in the form of this solitary wave was 
first obtained in [10] and later it was rediscovered a lot of times. 
Let us formulate the following theorem. 

Theorem 3.1. Let Y{z) he a solution of the equation 



y. + r^-^^o. (3.5) 

Then 
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+ i^y(z)-12/3y(zf (3.6) 

is a solution of the equation (|3.3|) in the case, 

- T - 625^ (^-^^ 

Proof. At the first step we should find the Newton polygon Li correspond- 
ing to equation (|3.3(l . For monomials of this equation wc have points: Mi — 
(-2,1), M2 = (-1,1), Ms = (0,2), M4 - (0,1), A/5 = (0,0). The support 
of equation H3.3|l is defined by four points: Qi — Mi, Q2 = M3, Q3 = M5 and 
Q4 = M5. Their convex hull is the triangle Li (Fig. 1). This triangle contains 
three vertexes r^°^ = Qj {j = 1,2,3) and three edges F^^^ = [Qi,Q2], T^^^ = 

[Q2,Q3],4^^ = [gi,Q3]. 




Figure 1: Polygon corresponding to the differential equation H3.3|l . 

Solutions of equation (|3.3|l have the second - order singularity. Let us make 
the following transformation 

w{z)=An+AiY{z) + A2Y{zf, (3.8) 

where Y{z) satisfies the first - order differential equation and has the first - 
order singularity. 
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Substituting 13. 8|) into the equation (|3.3|) . we obtain 



M2[w[Y]] = {(5Ai+2[3A2Y)Y,, - [u A^ +2u A2Y)Y,+ 




(3.9) 



The following points correspond to the monomials of this equation: Mi — 
(-2,1), M2 = (-2,2), Mg = (-2,2), M4 = (-1,1), Mr, = (-1,2), Mg = (0,0), 
My = (0,2), Afg - (0,4), Afg = (0,1), Mio = (0,2), Afn = (0,3), M12 = (0,0), 
A//13 = (0,1), Afi4 = (0,2), Afi5 = (0,0). The support of equation (EH is 
determined by seven points: Qi ~ Mi, Q2 = M2 = M^, Q3 = Afg, Q4 = 
Mil, Q5 = My = Mi4 and Qg = Mg — M13. Their convex hull is the quad- 
rangle (Fig. 2). This quadrangle has four vertexes F^"^ = Qj {j = 1, 2, 3, 4) and 

four edges F^'^ = [Qi,Q2], F^'^ = [Q2,Q3], T^'^ = [Qs^Qj], F^^ = [Qi^Qr]- So 
we have constructed the Newton polygon ^2- 

Following our method we should find the Newton polygon L3 with some part 
of edges parallel to those of L2- Besides that when the polygon L3 moves along 
the plane, his vertexes should cover the support of equation (I3.9|l . 



Figure 2: Polygons of the differential equations H3.9|l and 13.10(1 . 

We see that the triangle with vertexes M4 Q5 and Q7 satisfy these re- 
quirements. Thus the simplest equation is the Riccati equation with constant 
coefficients 




-2 -1 



1 



Ei[Y] ^ Y, + Y^ -b = 0. 



(3.10) 



Substituting 



Y,{z) ^ Ei[Y]-Y^ + b 



(3.11) 
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into equation (|3.9|) and equating coefficients at powers of Y{z) to zero, yields al- 
gebraic equations for parameters A2, Ai, Aq, b and Ci. Solving these equations, 
we get 



^' 5 25/3' 100/32' 

(3.12) 

" 2 625/32 

and the relation 

M2[w{Y)] =REi[Y], (3.13) 
where i? is a differential operator. This completes the proof. □ 

From equality H3.13|l we see that if Y{z) is a solution of equation (|3.11|) . then 
w{z) in formula (j3.8(l is a solution of equation (|3.3|l . 

Hence we have found the solitary wave in the form of a kink [10] 

6i/2 3z/2 / { v{x-Cot + ipfi)^^^ 



^^) = Co + — -— 1 + tanh ± '^^'^ ^u..,v-u, ^ (314) 
^ ^ " 25^ 25/3 V I 10/3 ' 

Here (^o is an arbitrary constant. 



4 Exact solutions of the nonlinear fourth — order 
evolution equation 

Let us look for exact solutions of the following fourth - order evolution equation 
[33] 

Ut-2Ux Uxx - Uxx -2uul+ Uxxxx = 0. (4.1) 

Using the travelling wave reduction H3.2() and integrating with respect to z, we 
get 

w^zz - wl - w'^ Wz - Cqw + C'l = 0. (4.2) 

Later let us prove the following theorem. 
Theorem 4.1. Let Y{z) be a solution of the equation 

Yzz ~YYz- Co ^0. (4.3) 

Then w{z) — Y{z) is a solution of the equation H4.2|l . provided that Ci = 0. 

Proof. Let us find the Newton polygon that corresponds to the equation l|4.2|l . 
The following points are assigned to the monomials of this equation: Afi = 
(-3,1), M2 = (-2,2), M3 = (-1,3), Af4 = (0,1), M5 = (0,0). The carrier 
of the equation contains five points Qi — Mi, Q2 — A/2, Q3 = A/a, Q4 = 
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Mi Q5 = M^. Their convex hull is the quadrangle Li with four vertexes 



rf - g, (j = 1,2,3,4) and four edges T]"' = [Qi,Q3], T^'^ = [^3,^4], T^'^ = 

[(34, Qs], ri^^ = [QiiQs] (see Fig. 3). Solution of equation 14.2|l have the first 
order pole. Suitable polygon L3 we find without making the transformation. It 
is the triangle in Figure 3. Thus we set 



(1) 



.(1) 



.(1) 



wiz) = Y{z), 

where Y{z) is a solution of the second order equation 
E2[Y] = - aVY, - 6 = 0. 



(4.4) 



(4.5) 




-3 -2 -1 

Figure 3: Polygons of the differential equations (|4.2() and 14.5(1 . 

Substituting 14.4(1 and 

Y,, = E2[Y] + aYY, + b (4.6) 

into equation 1(4. 2() and equating coefficients at powers of Y{z) to zero yields 
algebraic equations for parameters a, b and Ci. Hence we get 

a = 1, 6 = Co, Ci = 0. (4.7) 

We also obtain the relation 

M3[wiY)] ^RE2[Yl (4.8) 
where i? is a differential operator. This completes the proof. □ 
Integrating equation 1(4.3(1 with respect to z yields the Riccati equation 
Co 



Y,- Y' 



z + 2 C3 = 0, 



where C3 is a constant of integration. Setting 
2*. 



Y{z) 



(4.9) 



(4.10) 



8 



in H4.9|l we get 

+ ^ - ) * = 0. (4.11) 



^ 2 

This equation is equivalent to the Airy equation. Its general solution is 

^{z) = C\ Ai{- 2-1/3 Cq/^ z + C4} + C2 Bi{- 2-1/3 ^1/3 ^ _^ ^^1^ ^4^^2) 

where Ai{(^) and Si(C) are the Airy functions and Ci, C2 and C4 are arbitrary 
constants. 

5 Exact solutions of the generalized Kuramoto 
— Sivashinsky equation 

Let us look for exact solutions of the generalized Kuramoto - Sivashinsky equa- 
tion, which can be written as 

Ut + au"^Ux + S Uxx + P Uxxx + 7 Uxxxx = 0. (5.1) 

Equation H5.1|l at m = 1 is the famous Kuramoto - Sivashinsky equation 
[34,35], which describes turbulent processes. Its exact solutions at (3 — were 
first found in [34] . Its solitary waves at /3 ^ were obtained in [10] and periodical 
solutions of this equation also at /3 ^ were first presented in [13]. Recently it 
was shown [36, 37] that this equation did not have other solutions except found 
before. 

Using the variables 



X = X\ t =t—, u =u, a = a = — , (5.2) 

V 7 7 Vjo 

we get the equation in the form (the primes are omitted) 



Ut + au"" Ux + Uxx + CrUxxx + Uxxxx = 0. (5.3) 

Using the travelling wave reduction 
u{x,t)=y{z), z — x — Cot (5-4) 
and integrating with respect to z, we get the equation 

MM = y,,, +<yv.. + - Coy + 2/™+i = 0. (5.5) 

m + 1 

Here a constant of integration is equated to zero. 
Let us present our result in the following theorem. 
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Theorem 5.1. Let y{z) be a solution oj the equation 



2to3 _|_ 11 ^2 _|_ 28„^-^_ 9 



y ^ T 



V2m2 + 18m + 27 



(5.6) 



Theny{z) is also a solution of the equation provided thatm ^ 0, m —1, 



m ^ -|, m 7^ -3, m ^ and 



3(2to2 + 9to + 9) 



cr = ±- 



3(3 + m) 



V2m2 + 18m + 27 



(5.7) 



(2to2 + 18 77i + 27)^/^' 

Proof. Monomials of equation H5.5() are determined by the following points: 
Ml = (~3, 1), Ah - (-2, 1), Ma = (-1, 1), A/4 = (0, 1) and M5 = (0,m + 1). 
The Newton polygon Li that corresponds to equation H5.5|l is the major triangle 
in Figure 4 ((a): m > and (b): m < 0). Suitable polygon L3 can be found 
without making the transformation. It is the smaller triangle in Figure 4. This 
triangle containes the points: Qi ~ (—1,1), Q2 — M4 = (0,1), and Q3 — 
(0,1 + f). 




1 +m 



(a) 



(b) 



Figure 4: Polygons of the differential equations ()5.5|) and H5.8|l at m > (a) 
and m < (b). 

Consequently the simplest equation corresponding to L3 can be written as 

Ei[y] ^ y.-Ay"^ -By = 0, (5.8) 

where A and B are parameters to be found. 
The general solution of H5.8|l takes the form 



viz) 



Ci exp 



B m z 



(5.9) 
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Substituting 

= EM+ Ay"^ + By (5.10) 

into equation H5.5|l and equating coefficients at powers of y{z) to zero yields 
algebraic equations for parameters A, B, Cq and a in the form 

(2 m + 3) (m + 3) (to + 1) + 9 a = 0, (5.11) 

{m + l){B^ + aB^ + B - Co) ^ 0, (5.12) 

A^ (3 + m) (to + 1) {Bm + 3B + a)^ 0, (5.13) 

(to + 1) (2 cr^TO^ + 18 CT^ TO + 27 - 9 - 54 TO - 81) = 0. (5.14) 
Equations H5.11|l . (|5.12() . H5.13|l . and (|5.14|l can be solved and we obtain 
A = {-9ay/^ {m^ + 11771^ + 18m + 9y^^^ (5.15) 
(In final expressions we will take into consideration Ai only.) 



Co ^ B (B^ +aB + l) , B=-- , (5.16) 

^ ' 3 + TO 



3(3 + TO^_ (5.17) 
V2m2 + 18TO + 27 

So we have found parameters of the equation (|5.9|) , conditions H5.7|l and the 
relation 

M:,[y]=REM, (5-18) 

where i? is a differential operator. This completes the proof. □ 
The general solution of equation H5.7|l can be presented in the form 



y(z) = (C + Ci expi±— "^^ =i) (5.19) 
' V I V2TO2 + 18TO + 27jy ^ ' 

where Ci is an arbitrary constant and C is determined by expression 



C=± J0^^^2^!±I!^±^. (5.20) 
3^(to + 3)(2to2 + 5to + 3) 
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We would like to note that the value of solitary wave velocity H5.19|l tends 
to Co = T l/"\/27 as m 0, but at the same time Co ^ as to ^ oo. 
Assuming to = 1 in (|5.19() we obtain the known solitary wave 



and the value of the parameter a: a — ±12/v47. This solution of the Ku- 
ramoto - Sivashinsky equation is the kink [10-14]. 

6 Self similar solutions of the fifth — order Ko- 
rteveg — de Vries equation 

Let us find exact solutions of the fifth - order Korteveg - de Vries equation 
using our approach. This equation can be written as 

Ut + U^xxxx - 10 U Uxxx - 20 UxUxx + 30 = 0. (6.1) 

It has the self similar solution 

1 X 

where y{z) satisfies the nonlinear ODE of the form 

Vzzzzz - 10 y Hzzz - 20 + 30 - z - 2 ?/ = (6.3) 

Our results are summarized in the following theorem. 
Theorem 6.1. Let Y{z) be a solution of the equation 

Yzzzzz - 40 YY,, Y, - 10 y2 Y^^^ _ lo + 30 Y^Y, ~Y~zY,^0 (6.4) 
Then 

y{z)=±Y,-Y^ (6.5) 
is a solution of equation (|6.3|l . 

Proof. The following points correspond to the monomials of this equation: Mi — 
(-5,1), M2 = (-3,2), M3 = (-3,2), M4 = (-1,3), M5 = (0,1) and Mq = 
(0, 1). These points generate the support of equation H6.3|) : Qi = Mi = (—5, 1), 
Q2 = Mi = (-1,3), Q3 = M5 = Me = (0, 1) and Q4 = M2 = M3 = (-3,2). 
Now we can plot the Newton polygon Li (see Fig. 5) 

Solutions of equation (|6.3(l have the second order pole. Thus we can make 
the following transformation 

y = AiY + A2Y^ + A3Y,, (6.6) 
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Figure 5: Polygon corresponding to the differential equation of ijOJ- 

where Y{z) is a function of the first order pole. Substituting 1)6. 6|l into the 
equation studied, we obtain 

-^lYzzzzz + ^2^^^^^^ + -^'iYzzzzzz^ 

-10 (AiF + AaF^ + ^3^^.) [A^Yzzz + + ^3^,...) - 

-20 {A^Yz + + ^3^.) (Ain, + A2YI + ^3^..) + (6.7) 

+30 {A^Y + ^2^2 + A:iYzf {A^Yz + ^2^' + A^Yzz) - 



- A^Yzz 


) -2AiY- 




2A3y. = 








Qe «2' 
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r r r 


/ /a 
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Qi 
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1 1 
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Figure 6: Polygons of the differential equations (|6.7|) and H().8|) . 

The support of this equation is defined by ten points: Q\ = (—6, 1), Q2 = 
(-5,2), Q3 = (-4,3), Qi = (-3,4), Q, = (-2,5), Qg = (-1,6), Qj = (0,2), 
Qs = (0,1), Qg = (—1,1), Qio = (—5,1). In this case, the polygon L2 is 
the quadrangle presented in Figure 6. Studying this polygon we can find the 
polygon L3 for the simplest equation. It is the triangle also presented in Figure 
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6. Consequently the simplest equation for Ht).7|l is 

E5[Y] = r^zzzz + mi YY,,Y^ + Y^Y,,, + mg Y^^+ 

(6.8) 

+TO4 Y^Y,^ +m5Y + mGzY^ ^0, 

where mi, 7712, ma m4, ms and mg arc unknown parameters to be found. 
Substituting 

Fzzzzz = £^5 [i"] - mi YY^^Y, - m2 Y^Y,,, - mg F/- 

(6.9) 

— m4 Y'^Yz — m5 y — mg z Y^ 

into equation H6.7(l and equating coefhcients at different powers of Y(z) to zero 
yields algebraic equations for parameters mi, TO2, ma, m4, and mg. As a 
result we get 

mi = —40, m2 = —10, mg = —10, m4 = 30, 

(6.10) 

ms = — 1, mg = — 1, 
and the relation 

M5[y{Y)] =RE5[Y], (6.11) 

where is a differential operator. This completes the proof. □ 

Equation H6.4|l can be integrated in z and we get the fourth - order analogue 
to the second Painleve equation 

Yzzzz - 10 Y^Yzz -lOYYz^ + 6Y^ - zY - /3 = 0, (6.12) 

where (3 is an arbitrary constant. Thus we have expressed solutions of equation 
()6.3(l (and consequently of equation H6.1|l 1 through solutions of H6.12(l . 

7 Self similar solutions of the fifth — order mod- 
ified Korteveg — de Vries equation 

Let us find self similar solutions of the fifth - order modified Korteveg - de Vries 
equation 

Ut - lOu^Uxxx - 'iOUxUxx - + SOu'^Ux + Uxxxxx = 0. (7.1) 

It has the self similar solution [38-40] 

u{x,t) ^ {5t)-^/'^w{z), z = x(5t)~l/^ (7.2) 
where w{z) satisfies 

M^w] = Wzzzz - 10 w'^Wzz - 10 wwl + 6w^-zw-/3 = 0. (7.3) 
Let us present our result in the theorem. 
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Theorem 7.1. Let Y[z) he a solution of the equation 



y,,, + 26yn,-6i;^-66^y^y, -363y4 + ^ = 0. (7.4) 

Then 



w{z) = ±hY{z) (7.5) 
is a solution of the equation (|7.3I) . provided that /3 ~ Tl/2. 

Proof. To begin with the Newton polygon corresponding to the equation (|7.3() 
should be found. The following points: Ah = (-4,1), Ah = (-2,3), M3 = 
(-2,3), M4 = (0,5), M5 = (1,1), Mq = (0,0) are assigned to the monomials 
of the studied equation. The support of the equation contains five points Qi — 
Ah, Q2 = Ah, Q3 = M5, Qi = Me Q5 = Ah = Ah- Their convex hull is 
the quadrangle Li with four vertexes F^^'' = Qj (j = 1,2,3,4) and four edges 

rl'^ = [Qi,Q2], 4'' = [Q2,Q3], 4'^ = [Q3,Q4], ri^^ = [gi,Q4] (fig. 7). 
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r jfi rv 
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-\ ^3 - 


2 -1 





Figure 7: Polygons of the differential equations (|7.3|) and H7.8|) . 

Leadings members of the equation (|7.1|l are located on the edge F^^'' = 
[Qi, Q2]- Substituting w{z) = oq z^ into the equation 

w,,,,-\Ow'^w,,~lOwwl + Qw^ ^Q, (7.6) 

we get fours families of power asymptotics for solutions of H7.3|l : (p, oq) = 
(— 1, 1), (—1, —2), (—1, —3), and (—1, 4). Again taking into account the pole 
order of solutions and the results of section 2, we can look for exact solutions 
of the equation H7.3|l in the form 

w{z)^AiY{z), Ai^O, (7.7) 
where Y{z) is a solution of the following simplest equation 

Y,,, -aYY,, -bY^- cY^Y, -dY^-mz = 0. (7.8) 
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It is important to mention that the transformation (|7.3|) does not change the 
quadrangle Li and the Newton polygon L3 for the simplest equation H7.7|l pos- 
sesses all necessary properties. In other words Li = L2 and L3 is the triangle 
(see Fig. 7). The parameters a, b, c, d, and m can be found. 
Substituting H7.7(l and 

Yzzz = E3[Y] +aYYzz + bY^ + cY^Y, +dY^ + mz^O (7.9) 

into equation l|7.3|) and equating coefficients at different powers of Y{z) to zero 
yields algebraic equations for parameters a, b, c, d, m, /3 and Ai. Finally, we 
obtain 



-2 6, c = 66^ d = 3b^, to = — \, 

2b 



A,,2 = ±b, /3 = T^ 



(7.10) 



and the relation [39] 

M4w{Y)] ^RE^iY], (7.11) 
where i? is a differential operator. This completes the proof. □ 

Consequently there exist solutions of the studied equation expressed through 
the function Y{z) at /? = 1/2 or /3 = —1/2. Let us show that solutions of the 
equation H7.4|l are associated with solutions of the first Painleve equation 

Vzz^Cv^ +Dz. (7.12) 

Making the transformation 

v{z) = Aw^ + Bwz (7.13) 
in the equation H7.5|l . yields the equation 

Bwzzz + 'iAwwzz + {2A-CB'^)wl-2C ABw'^Wz-C A^w^-Dz = 0. (7.14) 

We see that the Newton polygon corresponding to this equation coincides with 
the triangle L3 at Figure 6. Combining H7.4|l and H7.1()|l . we obtain 

w{z) = Y{z), (7.15) 

provided that A = 6, i? = 1, C = 36, and D = —1/26. In other words, if Y{z) 
is a solution of (|7.4|) . then v{z) — bY'^ + Yz is a solution of 

Vzz^ibv'' - ^z. (7.16) 
We would like to mention that without loss of generality it can be set 6 = ±1. 
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8 Solitary waves of the sixth order nonUnear 
evolution equation 



Let us find exact solutions of the sixth order nonlinear evolution equation that 
is applied at description of turbulent processes [41,42] 



Ut + UUx + fi Uxx + <5 Uj^xxx- + £ Uxxxxxx = 0. (8.1) 

It is known that Kuramoto - Sivashinskiy equation and the Ginzburg - 
Landau equation that are used at description of turbulence are non-integrable 
ones because they do not pass the Painleve test [5]. However, these equations 
have a list of special solutions [5, 10, 13,36,38]. 

Equation 1)8. l|l does not pass the Painleve test as well and thus this equation 
is also non-integrable. However one can expect that equation (|8.1() has some 
special solutions. This equation admits the travelling waves reduction 

u{x,t)=y{z), z ^ X — C'ot, (8.2) 

where y{z) satisfies the equation 

Ci - Coy + ^j/^ + f3yz + Sy^zz + eyzzzzz = 0. (8.3) 

Ci is a constant of integration. As a result we get the following theorem. 
Theorem 8.1. Let Y{z) be a solution of the equation 



Yz+Y^-ak=0, (fc = 1,2,...,6). (8.4) 
Then 



/2520 

y{z) = 30240er^ + —— (5 - 50400 e ) 



V 11 



2520 , ^^^^^ 2 1260 ^ 12600 ^ 
^^.. + 20160s..^ + ^/?-^-jy + Co 

is a solution of the equation (|8.3() provided that 

^ 4112640 js^WfeS 5080320 J^^fe^^ 
C, = ~ 9999360 — ^ 



(8.5) 



55460160 5^wk^ 660240 S'^Wk^P 25200 5^Wk^ 
30371 ~^ 2761 ^2 — 5 



I2m P'^Swk , 12600 135^ Wk ,1^2 



251 



213811840 e^a^s _ ;L0204656(5£2afe2 _ 2045 (5^ - 92400 (J^eafe 
^~ 12l£ (9240eafc + 79^) ' 
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where ak and Wk are found from 



ak = (fc = 1,2, ...,6), 



W-l = , W2 = , W-t = , 

220 176' 440' 



(8.9) 



1 46031 
= 52800 '^'^ 



(8.10) 



1 / 46031 m ^ i%/3 / 46031 \ 

W5 6 = h 557 ± m H 

' 52800 y 2m 2 2 \ m J 

m = (113816753 + 1260%/8221079733)3 « 610, 966. 



(8.11) 



.12) 



Proof. The following points correspond to the monomials of equation 1)8. 3|l : 
Qi = (-5,1), O2 = (0,2), Q3 = (0,0), Q4 = (-3,1), andQs = (-1,1)- The 
Newton polygon Li of the equation is the triangle (Fig. 8). 
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Figure 8: Polygon corresponding to the differential equation H8.3|l . 

Solutions of equation (|8.3|l have the fifth ~ order pole. Let us express them 
through a function Y{z) possessing the first - order pole 



y{z) = Ao + AiY + A2Y^ + A^Y^ + A^Y^ + AgF^ 



.13) 



Substituting transformation H8.5() into (|8.3|) . we get new equation and the cor- 
responding polygon L2 in the form of quadrangle (Fig. 9). One of the suitable 
polygons L3 is the triangle in Fig. 9. This triangle is assigned to the Riccati 
equation with constant coefficients 



Ei[Y] = Y^ + Y^ ~a^0. 
Substituting 



(8.14) 



.15) 
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Figure 9: Polygon corresponding to the transformed differential equation H8.3|l 
and to the simplest differential equation (|8.14|l . 

into equation H8.4|) and equating coefRcients at different powers of Y{z) to zero 
yields algebraic equations for parameters A5, ^4, A3, A2, Ai, Aq, Ci, /3, ak and 
Wk- Solving these equations we have 



25205 

30240 £, ^4 = 0, A3 = — 50400£a, ^2=0, 



Ai 



2520 
11 



(5 a + 20160 £ 



11 
1260 



(8.16) 



251 



^ 12600 , _ „ 

^" 30371 7' 



For the parameters Ci, /3, ak, and Wk we obtain expressions (|8.6|l . (|8.7|l . 1)8. 8|l . 
(18.9(1 . H8.10|l . and (|8.12(l . We also see that the following relation 

M5[yiY)] ^RE,[Y], (8.17) 

holds, where i? is a differential operator. This completes the proof. □ 

The general solution of equation 1(8. 4|l is the following 

Y{z)^ ^tSinh{^z + ipo), (fc = 1,2,...,6). (8.18) 

Substituting ((8.181) into ((8.5(1 and taking into account that ak — WkS/e 
{k — 1, 6) we get six different solitary waves of equation ((8.3(1 . 
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9 Exact periodic solutions of equation ( 18.11 ) 

It was mentioned in the previous section that solutions of equation 1)8. 3|l have 
the fifth order pole. Thus let us make the following transformation 

y (z) = Bi+ B2R{z) + Bgi?^ + B^R^ + B-^RR,, (9.1) 

where Bk (fc = 1, 5) are constants and R = R{z) is a function of the second 
order pole. Now let us prove the following theorem. 

Theorem 9.1. Let R{z) be a solution of the equation 



2 1 RS^ , 1 RV21S^ 



Ri ^~-2R-' + aR-' - -a^R H — ± — + 

6 726 £2 2541 



1 3 , 13 , 1 V216^ 1 aS^ _ 1 aV2lS^ 



(9.2) 



108 " ^ 359370 ^ 119790 4356 ^ 15246 

Then 



y(z) = Co + 630 ( ea + — - 6eR ] R^ (9.3) 



11 

solution of equation iH.l\i .vrovided that 



(1,2) _ 10854 , 1^2^ 2484 V2lS^ 
^1 "161051 + 2^° ^ 161051 ^ ' 

Proof. Substituting 1)9.1(1 into equation (|8.3(l we get a new equation of the sixth 
order. The Newton polygon L2 of this equation is presented in Figure 10. Again 
we should construct suitable polygon L3. Let us take the triangle corresponding 
to the equation for the elliptic function 



El [R] ^ Rl + 2i?3 - ai?2 -2bR-d = 0. (9.6) 

This triangle satisfies necessary requirements. Substituting 
Rl ^ Ei[R] - 2R^ + aR^ + 2bR + d, (9.7) 



R.. - ^ -3R' + aR+b (9.8) 

2Rz az 
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Figure 10: Polygons corresponding to the transformed differential equation (|8.3fl 
and to the simplest differential equation (|9.t)|l . 

into equation (|8.3() and equating coefhcients at different powers of R(z) and 
to zero yields algebraic equations for the coefhcients i?5, ^4, Ba, i?2, Bi, and 
_Bo. As a result we obtain the following expressions 



B5 = -3780e, B4 = 0, B3 630£ a + ^ 5, 



-62 — 0, -Bi — Co 



(9.9) 



^ 121 e ' 



(9.10) 



1 2 1 52 1 

Oi 2 = a H ^ ± T. — , 

12 1452 £2 5082 £2 ' 



(9.11) 



1 3 13 J'^ 1 V2T(53 
^'^^lOs" 359370 £3 119790 ~ 



1 a(52 1 a\/2T(52 



4356 £2 ' 15246 £2 
and the relation 



(9.12) 



M^{y[R)\ =REi[R], 

where i? is a differential operator. This completes the proof. 



(9.13) 
□ 
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If i?2, and i?3 such that Ri > R2 > R3 are real roots of the equations 



2R^-aR^+r-a^-^^T^AR-^a^- 
\ 6 726e2 ^ 2541^2 I 108 



13 53 ^^2 a(52V2T „ 



359379 £3 119790e3 4356e2 15246£2 
then the general solution of 19.2(1 can be written as 



Ri^)^R, + iR,-R,)cnHzJ^^,Sl = (9-15) 

Thus we have found several solutions of equation 1)8. l|l at different values of 
the parameters. These solution are solitary and periodic waves (see expressions 
JHSl and (FOll ^. 



(9.14) 



10 Conclusion 

In this paper a new method for finding exact solutions of nonlinear differential 
equations is presented. Our aim was to express solutions of the equation studied 
through solutions of the simplest equations. To search the simplest equation we 
applied the power geometry and constructed the Newton polygons of nonlinear 
differential equations. The method is more powerful than many other methods 
because the structure of a solution is not fixed but can be found using the 
stated algorithm. Thus one can obtain quite general classes of exact solutions. 
It is important to mention that transformations between nonlinear differential 
equations can be also found with a help of our approach. As an example of 
our method application exact solutions of some nonlinear differential equations 
were found. In particular, solutions of the fourth - order evolution equation 
(14.11) were expressed in terms of the Airy functions; the connection between self 
- similar solutions of the fifth - order Korteveg - de Vries equation (|6.1|l and the 
fifth - order modified Korteveg - de Vries equation 1(7.11) were constructed; one - 
parametric family of exact solutions for the generalized Kuramoto - Sivashinsky 
equation 1(5.1(1 was found. Besides that we obtained solitary waves and periodic 
solutions of the sixth - order evolution equation ((8.1() that is used at description 
of turbulence [42,43]. 
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